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ABSTRACT 

Noncommutative geometry can provide effective description of plrysics at very sliort distances talc- 
ing into account generic effects of quantum gravity. Inflation amplifies tiny quantum fluctuations in 
the early universe to macroscopic scales and may thus imprint high energy physics signatures in the 
cosmological perturbations that could be detected in the CMB. It is shown here that this can give 
rise to parity-violating modulations of the primordial spectrum and odd non-Gaussian signatures. 
The breaking of rotational invariance of the CMB provides constraints on the scale of noncommuta- 
tivity that are competitive with the existing noncosmological bounds, and could explain the curious 
hemispherical asymmetry that has been claimed to be observed in the sky. This introduces also 
non-Gaussianity with peculiar shape- and scale-dependence, which in principle allows an independent 
cross-check of the presence of noncommutativity at inflation. 

Subject headings: Cosmological inflation, Statistical anisotropy of the cosmic microwave background, 
Scale-dependent non-Gaussianity, Noncommutative spacetime 



1. INTRODUCTION 

The statistics of the temperature anisotropics in the 
cosmic microwave background (CMB) are measured by 
the Planck satellite to an unprecedented accuracy. This 
allows to efficiently probe, in addition to higher-order 
correlations, i.e. possible non-Gaussianity, the detailed 
structure of the two-point correlations, i.e. possible sta- 
tistical anisotropy. 

There is a number of anomalies already in the present 
data, which have raised a lot curiosity both from the 
theoretical side as well as from the data analysis side 
(Hansen et al. 2008; Hoftuft et al. 2009; Koivisto & Mota 
2006, 2008a,c; Rakic et al. 2006; Koivisto et al. 2009, 
2010; Mota et al. 2007; Rakic & Schwarz 2007; Koivisto 
& Mota 2008b; Jimenez et al. 2009; Zumalacarregui et al. 
2010; Copi et al. 2010; Bennett et al. 2010). In particular, 
the hemispherical asymmetry, first reported by Eriksen 
et al. (2004), seems a quite unexpected feature within 
the standard model of cosmology and hasn't yet been a 
satisfactorily traced to a possible systematic error. The 
question whether the universe is odd was asked in (Land 
& Magueijo 2005), and there are recent investigations 
(Kim & Naselsky 2010a, b) finding hints of evidence for 
a positive answer. 

This prompts to look for possible cosmological ori- 
gins of odd-parity statistical anisotropics. In the present 
study, we investigate the effects of noncommutative ge- 
ometry to the primordial spectrum of perturbations, usu- 
ally assumed to be generated by quantum effects during 
infiation or shortly afterwards. The observational impli- 
cations are derived, in terms of the harmonic coefficients 
of the CMB spectrum, and the non-Gaussianity param- 
eter /tvl- We find that in general the noncommutativ- 
ity of spacetime geometry induces parity violating mod- 
ulations of the spectra of fluctuations, thus generating 
distinct signatures in the statistics of CMB. In particu- 
lar, this suggests that the hemispherical asymmetry (and 
various other anomalies) could originate from the funda- 



mental properties of spacetime that are relevant at the 
vast energy scales at play in the inflationary epoch. 

Chu et al. (2001) remarked that spacetime noncom- 
mutativity can be constrained by the statistics of in- 
flationary fluctuations. The power spectrum has been 
computed (Tsujikawa et al. 2003; Koh & Brandenberger 
2007) and various other aspects of noncommutative in- 
flation have been discussed in the literature (Alexander 
et al. 2003; Calcagni 2004; Palma & Patil 2009). The 
CMB constraints beyond the power spectrum have been 
explored also (Lizzi et al. 2002; Akofor et al. 2008a, 2009; 
Karwan 2010). Here will adopt the formalism of Akofor 
et al. (2008a). What is new in particular, is that we 
point out the presence of odd signatures and compute 
the structure of the two-point and three-point correla- 
tors in more detail and generality than previously. 

We also clarify an ambiguity of the results, which forces 
us to introduce an additional parameter. To assess the 
robustness of the results, we consider in addition the al- 
ternative approach of Kobakhidze (2008). It becomes 
clear that the details of the predictions can depend upon 
the particular model, but there are generic features which 
appear already in the simplest cases (in particular, in the 
case of canonical noncommutativity with constant 9 in 
the comoving frame). These nontrivial statistical fea- 
tures may thus be present, at an observable level, even 
in the simplest inflationary ACDM models, if one takes 
into account the effect of spacetime uncertainty principle 
on the inflationary fluctuations. 

However, other means of generating parity violations 
can be introduced too. A simple way is to assume an 
inhomogeneity present at the early universe. (Erickcek 
et al. 2009) considered that a large scale perturbation 
of the curvaton field might result in a power asymme- 
try. Tangen (2009) calculated the CMB pattern from 
a single superhorizon perturbation, which indeed shows 
couplings between adjacent multipoles. This is different 
from the approach of considering dipole in the primor- 
dial spectrum, which introduces adjacent-mode correla- 
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tions for the anisotropies of the random fluctuations at 
all scales, as will become clear below. One may also 
contemplate on possible parity-violating couplings of the 
inflaton field. Alexander (2008) has considered the possi- 
ble role of Chern-Simons terms (Chern & Simons 1974). 
Finally, spontaneously broken isotropy, occurring due 
to imperfect dark energy, has been shown to produce 
odd modulations (Gordon et al. 2005; Axelsson et al. 
2010). There are qualitative differences to the present 
case, which will be clarified in section 5.3. 

In the following section 1 we review the basic results 
of inflationary perturbations and discuss how these can 
be applied when the spacetime is noncommutativc. We 
then implciment this in section 3 in the case of canoni- 
cally deformed spacetime commutation relations and in 
section 4 in a framework based on deformed Heisenberg 
algebra of quantum fields. We are then ready to discuss 
the observable patterns in the CMB sky. The properties 
of the two-point functions and of the non-Gaussianities 
are clarified in section 5, and section 6 is a brief con- 
clusion. The CMB two-point correlation in terms of the 
multipole expansion of the primordial spectrum is given 
in the section 5.1. 



commuting case. On superhorizon scales we can write 

Ce {t, x) = N'{t, ti)5(t)0 {ti, x) + \N"{t, ti)5<t>e {ti,xf + ... , 

(1) 

where the subscript 9 is introduced to denote noncommu- 
tativc variables. The scalar field perturbations 5(j)e gen- 
erated during inflation are evaluated on a uniform curva- 
ture hypcrsurface ti soon after the horizon crossing of all 
the modes of interest. Their properties differ from the 
corresponding commuting quantities 5(f)Q as we will dis- 
cuss below. The function N{t, U) measures the number 
of e-foldings of a classical Priedmann-Robertson- Walker 
(FRW) universe from the uniform curvature hypcrsur- 
face at ti to a uniform density hypersurface at some final 
time t when the universe is evolving adiabatically. The 
primes denote derivatives with respect to the classical 
background field <j>. The derivatives of N{t,ti) describe 
entirely classical properties of the theory and their val- 
ues coincide with the corresponding commutative theory. 
The curvature perturbations produced in a noncommu- 
tativc and commutative theory with the same classical 
solutions therefore differ only by the different properties 
of 5(1)0 and 5(j)o. We turn to discuss the relation between 
S<p0 and in more detail after briefly reviewing some 
standard results for S^q. 



2. CURVATURE PERTURBATION IN 
NON-COMMUTATIVE INFLATION 

In the vast majority of models, primordial perturba- 
tions originate from quantum fluctuations of light scalar 
fields produced by the inflationary expansion. Their 
properties depend on the physics operating at the very 
high energy scales present during inflation. It is conceiv- 
able that physics at such high energies becomes inher- 
ently non-local; such models can be effectively described 
by noncommutativc theories. Here we are interested in 
studying primordial perturbations generated in noncom- 
mutative theories of inflation. We treat gravity as a clas- 
sical background which is not affected by the noncom- 
mutativc effects. In this approach the noncommutativ- 
ity affects only the properties of quantum fluctuations 
generated during inflation. 

To keep the discussion transparent, we restrict our 
analysis on general single fleld models where primordial 
perturbations effectively arise from fluctuations of a sin- 
gle scalar degree of freedom while additional scalars 
may affect the background evolution. This class of mod- 
els obviously contains the standard single field infiation 
in which case (f> is the inflaton. In general, however, cp 
can be a scalar fleld different from the inflaton-like fields 
which dominate the energy density. Well known exam- 
ples are the curvaton model and modulated reheating 
scenario where primordial perturbations can arise solely 
from fluctuations of a light field (p which remains sub- 
dominant during inflation but affects the expansion his- 
tory at a later stage (Lyth et al. 2003; Dvah et al. 2004). 

The primordial perturbations are conveniently charac- 
terized by the curvature perturbation ( which measures 
fluctuations in the spatial curvature on uniform energy 
density hypersurfaces. Since we take gravity as a classical 
background which is not affected by the noncommiitative 
effects, the curvature perturbation can be computed us- 
ing the i5A'' formalism in close analogue to the standard 



2.1. The mode functions in the commutative case 

We consider the FRW metric in terms of conformal 
time T and including scalar perturbations in the Newto- 
nian gauge in the absence of shear (Mukhanov 2007): 

ds^ = -a2(T) [dr^ (1 + 2$) - da;2 (1 + 2$)] . (2) 

A scalar field ip can be expanded in terms of the annihi- 
lation and creation operators as 

i.{x) = I («.(r)e*-a.+t.Kr)e-*-4)-^^|^. 

(3) 

The operators satisfy the canonical commutation rela- 
tions. The canonical momentum can then be identified 
as 

7r(x) = ^ (a(r)V(x)) ^ {a^j)' , . (4) 
mode functions w(r) obey the time evolution equation 



Wfc + ( fc^ - — ) Mfc = , 



(5) 



with the well known Hankel function solutions that, when 
matched with the initial Bunch-Davies vacuum at early 
times, reduce at late times outside the horizon to 



g-ik(T-Ti) / j 

Ukir) = ^= — 7- - 1 

jl.p-ik{T-Ti) / i I 



2k 



kr k'^r'^ J 



(6) 

(7) 



The spectrum of scalar metric perturbation $ in the con- 
formal Newtonian gauge is then related to the spectrum 
of the field fluctuation as 



P^{k) = ^^P^{k)., 



(8) 
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and evaluated at the horizon crossing a{ri)H = k, where 
the Hubble rate is approximately constant when the slow- 
roll parameter 

H 



(9) 



is small. The spectrum of the metric perturbation $0 in 
standard single field inflation can be known to be given 

by 

3. CANONICAL NONCOMMUTATIVITY 

A canonical way of deforming the spacetime is to intro- 
duce the commutation relations for the coordinate oper- 
ators 

[xt,,x^]=i9fj,^, (11) 

where in the simplest case Of^v is an antisymmetric con- 
stant matrix of dimension length squared. It is well 
known that this is the exact low-energy limit of open 
string theory with a constant antisymmetric background 
field (Seiberg & Witten 1999). In general, a commu- 
tation relation of the form (11) induces the uncertainty 
relation for coordinates 



Ax" Ax" > -I 



(12) 



so that a spacetime point is heuristically replaced by 
a Planck cell. The ordinary coordinates may then be 
thought to be obtained by coarse-graining over scales 
smaller than the fundamental scale of order -v/j^. Thus, 
noncommutativc spacetime provides a framework that is 
compatible with generic features of quantum gravity like 
the uncertainty principle and nonlocality. 

The commutation relations (11) generally assumes 
more complicated form when expressed in alternative co- 
ordinate systems. It is thereby essential to specify in 
which frame this relation is taken to hold as written down 
above. In cosmology, a natural frame to consider is the 
comoving one. We call the physical scale 9f/^, whereas 
9^1, is the matrix corresponding to the coordinates of 
an observer, to whom this matrix then is a constant 
throughout the evolution of the universe. We perform 
the computations in the comoving frame, but in the end 
translate the result into the physical scale employing the 
relations 

9^o1=a{t)9oi{t), 9^^ = a^t)9oi{t) , (13) 

where the 9^,y here and in the following is evaluated in 
the comoving coordinates. 

Consistent statistics in noncommutativc spacetime 
Akofor et al. (2008b) require deformation of the quan- 
tum fields by the exponential operator defined by the 
following relation: 



tpeix) = ipo{x)exp 



(14) 



where the lower index 9 refers to the deformation, so tpo 
is the corresponding field in the commutative case. In the 
following we will be interested in the two-point correla- 
tion of the inflaton field in this setting. This implies that 



the vacuum expectation value of the two-point function 
can now be written as 

{0\^l{x)Mx')\0) = e-i'"'-''"'''-^' {0\Mx)Mx')\0} . 

(15) 

Writing this in terms of the Fourier image ?!>e(k, t) defined 
by 

^</)(k,t)e*-, (16) 



we obtain 



(17) 



-i (fe.^-a„ -a,e-feO (OI0S (k, i)</'o (k, t') |0)e'- 



kx 



where we have used the fact that in the usual case ^ = 
(only) the different wavemodes are uncorrelated, 

(O|(/.S(k,t)0o(k',t)|O) = (27r)3p^„(fc,t)(53(k-k'), (18) 

where P^{k,t) is the power spectrum. Let us call the 
time-space components of the noncommutativity the 
three- vector 9 as 

9°' = 9' . (19) 

By comparing the form (17) to the usual case we readily 
infer that 



(0|<^T(k,t)^e(k,t)|0) = 

lim^e-^^^-''(^*+^*') X (O|4(k,i)0o(k,O|O) 



(20) 



t->-t- 



{27rfP^Jk,t~-9-k 



Consider then the spectrum of metric perturbation (10) 
in near de Sitter space where H is approximately con- 
stant is given by 



-1 



-Ht 



Ha{t) 



(21) 



Using this we can immediately combine equations (20) 
and (10) to obtain the spectrum in noncommutativc ge- 
ometry. Evaluated at the horizon crossing, we have 



P3,,(k) =P$„(A;)e 



(22) 



Thus the spectrum will be direction-dependent. Further- 
more, it is not parity invariant. We see that the leading 
correction is a dipole with an amplitude Aim ~ |^| and 
blue-tilted spectral index ni^^ ~ 2. The next correction 
is the even-parity quadropole term, with an amplitude 

-42, m ^ ^I^P and spectral index ni.m ~ 3, and so on. 

However, wc have to choose the consistent parts of the 
correlator in order to obtain a physical result. As clear 
from the subsection 5.1, the odd multipole modulations 
should have imaginary coefficients, otherwise the result 
is not sensible as the real-space correlators and the CMB 
sky would not turn out real. This would be cured if we 
promote 9 into an imaginary parameter but this seems 
inconsistent with 11. Therefore we adopt the following 
prescription 



(...) a{...)M + a){...)A, 



(23) 
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where (. . .) denotes schematically some a correlator, 
{. ■ ■) M is its self- adjoint, and {■ ■ ■) a its anti-self-adjoint 
part. Then a G [0, 1] is a parameter which corresponds 
to some kind of phase. Akofor et al. (2008a) considered 
only the self-adjoint part of the correlators, which cor- 
responds to a = 1 in our parametrization. However, we 
do not know any physical rciason why the remaining part 
should not contribute to observed correlations. As we are 
unable to determine the value of a from first principles, it 
is left as a parameter to be determined by observations. 
We then introduce the notation 

exp^(a;) = a cosh (a;) -|- i{l — a) sinh {x) . (24) 

In this prescription, the result (22) becomes 



P$,(k) = P$„(fc)exp„(if^-k), 



(25) 



which is the main result of this section. 

The conceptual problem noted above appears in alter- 
native frameworks too and is thus not merely a possible 
inconsistency in of the particular formalism employed 
here. In particular, one could also start by expanding 
the action in the noncommutative parameter and derive 
the correlations of the perturbations from the ensuing 
equations of motion. One obtains also then imaginary 
results, see e.g. the four-point hmction calculated by 
Chu et al. (2001). It can be pointed out also that the 
problematics of observing correlations of noncommuting 
observables are independent of the nature of noncommu- 
tativity. Noncommutativity between space and time has 
not yet been put into theoretically rigorous footing but 
generically seems to imply unitarity violations (Gomis 
& Mehen 2000; Chaichian et al. 2001). Though in open 
string theory with a constant electric background, which 
is supposed to exhibit noncommutativity between space 
and time, these problems are absent (Seiberg et al. 2000), 
the noncommuting field theory with constant ^^i, can be 
recovered from string theory only in the case of magnetic 
background field which then corresponds to vanishing 6. 
The need for the prescription (23) would nevertheless 
reappear for higher order correlations even when ^ = 
as will be seen in section 5.4. 

4. DEFORMATION OF THE HEISENBERG 
ALGEBRA 

Violation of microcausality in the spirit of (stringy) 
uncertainty principle can also be described by imposing 
noncommutativity of quantum fields (Dubovsky et al. 
2008) (instead of the coordinate operators, as in the 
previous subsection). Consider the following equal-time 
commutation relations in expanding spacetime 



[</'(x,T),(?i(y,r)]^ 

[7r(x, r),7r(y,r)] : 
[0(x,T),7r(y,T)]: 



a2(r) 



/(x-y), 



<5(x-y). 



(26) 

(27) 
(28) 



For notational convenience we parametrize /x(t) = 
/xoa(T)~, where the constant where /io is the character- 
istic scale of microcausality violation with the dimension 
dim[0]-dim[/]/2. The scale-factor dependence is added 
because we want to consider also the case where the form 



of the commutator is constant in comoving coordinates, 
n — 2. The time-dependence of the effective parame- 
ter /X does not affect the computation and we return to 
the different choices in the analysis of the results. The 
derivation here follows closely Kobakhidze (2008), where 
n = 0. The difference to the usual case is now only the 
odd function /(r) appearing in the first commutator. It 
is useful to note that by defining the field il^{x) as^ 



(l){x) = i){x) 



2a(T) 



J /(x-z)7r(z,T)d3z, 



(29) 



we recover canonical commutation relations for the pair 
(V',7r). This observation allows us, analogously to the 
previous subsection, to relate correlations in the noncom- 
mutative case to the correlators in the standard case. 
In particular, one may check that if the field tp satis- 
fies the equations for standard inflaton we described in 
section 2.1, one can translate those conventional results 
into the noncommutative set-up (26-28) by employing 
the shift (29). The noncommutative inflaton field is then 
expanded is terms of the mode functions Uk introduced 
in section 2.1 as 

Hx) =/ (0 {Mr) - ^F(k)<(r)) e*-a. 



(k)nr(r)) e-^^-ai 
where F(k) is the Fourier image of the function /(x), 

F(k) = J /(y)e*''yd3y (31) 

Since /(r) is an odd function, we have also F{—'k) = 
F(k). The two-point correlation function follows then 
straightforwardly: 

d^kUu{T)U\{T)e^<''-y'> 
a2(T) (27r)3 

:/p;(k)e^>^-(--)^. (32) 



(0|<^(x,r)</.(y,r)|0)= j 



where Ui^{t) is defined as 
C/k(r)=4(r)- 



M^(r) 



F(fc)4(T) . 



(33) 



In the second line we have identified the power spec- 
trum of the non-selfcommuting field (p and denoted it by 
P^(k). It is easy to see that lest physical observables 

become imaginary, the Fourier image F(k) must be real, 
and thus we cannot express the first line of (32) as a 
square^. Now using the relation of the infiaton and met- 
ric perturbation spectra (8) and the late-time limit of the 

^ The fcLctor l/a(T) was missing in Eq.(6) of KobaJchidze (2008). 

2 Our result (30) differs from the Eq.(13) in Kobakhidze (2008) 
by the sign of the argument of the second F(k.), but our (32) would 
agree with the Eq. (14) in Kobakhidze (2008) if F(k) was imaginary. 
This indeed seems to have been the assumption in Kobakhidze 
(2008), which then however leads to imaginary temperature corre- 
lations, see Eq.(23) there. 
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mode function solutions (6,7) wc obtain 
SttG 



1 + fcV^ (1 + if,''kF(k)) 



-|^F2(k)(l-fcV + fcV) 



(34) 



The leading order contribution in the parameter is thus 
odd in parity. 

Let us then look at some specific forms of the noncom- 
mutativity. A simple assumption for the form of the; func- 
tion /(r) is a delta- function and that the commutator is 
constant in comoving coordinates. An odd combination 
is 

M2(r)/(r) = AnHfiy [S{v - v) - J(-r - v)] . (35) 

In this prescription, the commutator (26) gets contribu- 
tion from spacelike separations equal to v. Then (34) 
becomes, neglecting the decaying modes 

P^(k) = P$(fc) [l + k^{^^oaTf sin(k • v) 



(36) 



Remarkably, at the leading order the predicted modula- 
tion has the same form as in the previous case (22), if 

we identify the vectors v — HO. One of course obtains 
hyperbolic sinus instead of the ordinary one by promot- 
ing V to an imaginary vector parameter (and dropping 
the i in (35)). The second line in (36) contains modes 
growing outside the horizon. They can be eliminated by 
choosing n < —1 in /i(T) = /ioa"(T). If n < —1 both the 
odd and the even contributions are decaying (regardless 
of the form of /(r)). The spectral index of the leading 
modification in (36) is ng -I- 3, but this depends sensi- 
tively on the precise form of /(r). As an example, the 
form 

TT^ r 

f{r)=i— ^{r^ - Vz)a{va, - rx)(T{vy - ry)a{vz - r^) 

+ irz +Vz)(t{vx +rx)(7{vy +ry)a{vz +rz) , (37) 

where a{x) is the sign of x, results in F(k) = sin(k • 
■v) / {kxkyki) , which results in a strongly blue-tilted spec- 
tral index. Thus we may obtain similar correlations as 
in section 3 by choosing a suitable function /(r). 
For the three-point function one gets 



{(t>{K)cj>{k2)<P{h)) = 



(38) 



(27ra(r))° 

>(ki + k2 - k3)[/k, (r) + ^^(ki - k2 - k3)f/^^ (r) 

in terms of the functions ?7j(r) defined in (33). 

5. PATTERNS FROM NONCOMMUTATIVE 
INFLATION 

In this section we discuss some observational implica- 
tions of the results at more length and derive an explicit 
expression for the non-Gaussianity. 



5.1. Multipole expansion of the primordial power 
spectrum 

The temperature anisotropy field is conventionally ex- 
panded in terms of the spherical harmonics and on the 
other hand considered in the Fourier space 



°° r d^k 

e(x,e,ry)=^ ^ agmYem = / -^r-^e*-^7eke(k,e,7?) 

e=om=-e ■' ^'^'^> 

(39) 

where we have normalized the transfer function 0(k, e, 77) 
with respect to the initial amplitude of the primordial 
curvature perturbation TZk). It follows that the coeffi- 
cients aim are given by 



^7ekF/„(A;)eKfc). (40) 



fk 
27r2 

where we have introduced the transfer function which 
depends only on the magnitude of the wavevector, 

&i{k)= I je{kr{v))e{k,v)dr] , (41) 



since we assume the evolution to be isotropic. The possi- 
ble anisotropy appears in the primordial spectrum, which 
can be expanded also in spherical harmonics as (we refer 
to this expansion of the primordial spectra always with 
capital letters indices to avoid confusion with the expan- 
sion of the temperature anisotropies) 



(7^k7^^) = 5^(k-k')|fp(k) 



(42) 



, 27r2V47r 



M 



fc3 



LM 



L=l M= 



ko 



ULM—l 



In the first equality we have used the WMAP conven- 
tions, and in the second one employed the parametriza- 
tion of Armendariz-Picon & Pekowsky (2009). They use 
the pivot scale ko = 2 ■ 10~3/Mpc. For the time being, 
we allow independent spectral indices ulm for each mul- 
tipole L,M. Using the formula (40) and the primordial 
spectrum (42) , the correlation matrix can be written as 



27r2 



(43) 

L=l M=-L 



We have separated here the integrated contribution from 
perturbations of all difi^erent magnitudes. 



LM 



tL,1\ 

leei 



dk f k_ 

k \ko 



Qe{k)&eik). (44) 



They are then weighted by the geometrical factors ^ 
which happen to be proportional to the coefficients of 
the Gaunt series. We may write them in terms of the 
Wigners 3-functions as 
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&'m' = (-ir+V(2^+l)(2^' + l)(2i+l) 
when £, i' and L satisfy the triangle condition. For the dipole these become 





I' 




( i I' L\ 









\ -m m' M I 



^L=l,M=-l _ /^r 



m',m— 1 



{(: + in - + m) - ni + l)[C - in + 2) 

-11/ ^,^n , <)^',^+l\ 



2(2^- 1)(2£+1) 



2(2£ + l)(2^ + 3) 



^L = l,M=l _ /^r 



m' ,m+l 



{e-m-l){e-m) , {£ + m + l){£ + m + 2) 

,^-11/ ^/^/l . 0£',^+ll 



2(2^- l)(2f +1) 



2(2^+ l)(2£ + 3) 



/(£-m)(£ + TO) - /(£-m+l)(£ + m+l) 



(2f- l)(2f + 1) 



(2£ + l)(2£ + 3) 



(45) 

(46) 
(47) 
(48) 



All odd multipole coefficients in the spectrum are imag- 
inary A21C+1 M — —■A2K+i,M, and the even are real 
•^2K M — ■^2K,M for any K. The geometric coefficients 
are symmetric, ^jj^^/^/ = ^.t'm'-im (o^® can check this 
is the case for the dipole above). The angular correla- 
tions of course turn out to be symmetric (a^m'^^'m') 
{ae'm'Ci}m), though the primordial spectrum may not be, 
(7^(x)7l(x)) 7^ (7?.(x')7?.(x)). As shown above, this can 
be understood by the noncommutative quantum nature 
of the fields whose fluctuations are responsible for the 
perturbations. 

5.2. Anisotropic power spectrum 

As found in sections 3 and 4, the two-point function ac- 
quires typically exponential modulations from noncom- 
mutative geometry. For clarity, let us focus on an ex- 
ponential term e'^ '"° in the following. This be decom- 
posed using the Rayleigh formula and by expressing the 
Legendre polynomial Pe in terms of sum of products of 
spherical harmonics: 

00 

gik.ro ^ ^ ^€(2£ + i)j,(fcro)P,(fc • fo) = (49) 

00 e 

e=o m=-e 

Now comparing the spectra with the general form (42), 
and using the orthogonality of spherical harmonics to- 
gether with (49), we obtain the amplitudes of each mod- 
ulations: 

47ri^ii(fcro)Y^V(firo)^- (50) 

Clearly the scale dependence of these coefficients cannot 
be described by simple power-laws. Instead, the modu- 
lations will be oscillating along the fc-modcs. 

Let us first comment the power spectrum, we note 
that already the isotropic spectrum is modified with re- 
spect to the usual result Aoo{^^' = 0) = ^, because 
of the nontrivial fc-dependence encoded in the function 
joikra) — sin(fcro)/fcro. In principle the oscillatory be- 
havior of the modulation could result in "wiggles" there 
seem to appear in the observed spectrum. Such wig- 
gles have also been predicted from transplanckian physics 
(Martin & Brandenberger 2003) or from cyclic inflation 



(Biswas et al. 2010), and the data has been shown to be 
compatible with such features. 

Furthermore, there is an inflnite series of higher- 
multipole modulations which will introduce statistically 
anisotropic correlations. The reflects the nonlocality of 
the underlying model. In principle all types of modu- 
lations are present, meaning that every £-modc; is cou- 
pled to any other. Each wavemode of perturbations that 
contributes to the power spectrum is also relevant to 
the anisotropic couplings. In particular, as one expects 
from ultraviolet noncommutativity, small wavelengths 
contribute most to the modulations at all scales. Contri- 
bution from extremely small wavelengths would cancel 
out due to rapid oscillations. 

It is crucial to note however, that in practise the scales 
at which the modulations are strongest, do not con- 
tribute to the angular modes observed in the CMB un- 
less the length scale of noncommutativity is much above 
the Planck scale. The spherical Bessel functions je{z) 
have their highest peak at about z ^ L, so each L- 
modulation will be strongest at wavemodes correspond- 
ing to k ^ L/rQ. On the other hand, the contribution to 
the CMB anisotropy Qe{k) from inhomogeneous sources 
Q{k, t) is also dependent on the spherical Bessel function 

0^(fe) = j Mkr)e{k,T)dT. (51) 

The scales contributing most to the multipole i are 
k ~ £/t^, where the comoving distance to the last scat- 
tering surface is about r, « 14000 Mpc. Thus the modu- 
lation of the order L will affect maximally the correlators 
corresponding to the multipole £ when 

ro^lA-W^''(^j^Mp'^H0^^, (52) 

where Mp^ is the Planck length and /i^^ ^ 6~i is the 
noncommutative length scale. So, the dipole modula- 
tion at the multipoles £ ~ 1000 would be of order one if 
jj. ~ 10"^° Mp ~ 10~^°Aqcd at inflation. Now, to trans- 
late this into the physical energy scale of noncommutativ- 
ity observable today in laboratory we should recall the 
relation (13). If the reheating temperature of the uni- 
verse was close to the GUT scale ~ 10^^ GeV, the scale 
factor at the end of inflation was about anH ^ 10~^^ 
of its value normalized to unity at the present. This 
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gives us fi 



ph 



10 TeV. We 



to 



note also that at those muhipoles the cosmic variance 
is negligible and Planck can be expected to measure de- 
viations from statistical anisotropy at percent level or so. 
Furthermore, the total effect of the modulations does of 
course not come from the peak of the Bessel functions in 
(50) but is the cumulative contribution integrated from 
all scales. We may then expect several orders of magni- 
tude improvement to the above estimate of the maximal 
noncommutative scale ji that may be observed in the 
CMB. Conservative lower bounds from modifications to 
standard model of particle physics give ^ > few TeV 
(Hewett et al. 2001; Mocioiu et al. 2000; Carroll et al. 
2001). Thus, the tightest bounds may turn out to be 
cosmological. 

Full comparison with the data would require consid- 
erable technical difSculties, firstly because all the ob- 
served multipoles should be included in the analysis, and 
even higher k-modcs than usually corresponding to those 
would have to be taken into account. Moreover, since 
there occur couplings between arbitrarily separated £- 
modes, one cannot employ the previous techniques that 
have been developed to deal with sparse correlation ma- 
trices (with only the diagonal and some adjacent entries 
nonvanishing) . Finally, the distortion of the power spec- 
tra should be tested in conjunction with the effects of the 
anisotropic correlations. 

Therefore, and because both observations and theory 
suggest these effects should be small, let us then, instead 
of the full pattern, consider the power series expansion 



l + zkTo-i(kTo)2 + ... (53) 



Note that the expansion of the more general parametriza- 
tion (25) is essentially very similar. It is useful to sepa- 
rate the magnitude r of ro = rf, defining the unit direc- 
tion vector f decomposed as 



r±=T 



\/2r 



To = 



r = \ro\ 



(54) 



Then the nonvanishing contributions to the spectrum 
may be written as the following. The amplitudes are 



(55) 



Aoo = A, = 2^-^/ —korf-A 



Aio = 2i^korroA , ^i(+i) = 2i^korr+A , 

and the corresponding spectral indices are 

noo = ns , nim = 1 + . (56) 

In a companion paper we test the leading order dipole 
correction given by (55) with the data, which is found to 
slightly prefer the presence of the dipole (Groeneboom 
et al. 2010). 

5.3. Comparison with imperfect source models 

Let us remark on the difference to the imperfect dark 
energy model, where similar geometric modulations can 
appear as well. There the dipole is due to an anisotropic 
source, it's contribution with respect to the quadropole 
is subdominant. This is because there then occurs can- 
cellation to the odd correlations, as they are proportional 



{aema*i+L)m')DE 



Ao J d{\ogk) [Qi{k)@t+i^{k) - Qi+L{k)et{k)] ,(57) 

where A quantifies the magnitude of the anisotropy, 
Sf{k) is the anisotropic transfer function and L is odd. 
The function (51) gives the source contribution in the 
isotropic case, in the presence of imperfect sources the 
O's depend also on the direction of the wavevector. Due 
to the partial cancellation effect in such a case, Axelsson 
et al. (2010) found the quadropole contribution dominant 
even though it was suppressed by the small parameter 
corresponding to Aq. However, since in the present case 
the dipole is of primordial origin, we have 



Alm / d{\ogk)Q,{k)Qt+L{k) , (58) 



{l+L)m'l 



and the magnitude of the odd and even modulations is 
expected to be similar. A priori, higher multipoles are 
again suppressed by some small parameter, and the dom- 
inating correction to the monopole is now generically the 
dipole. 

5.4. Non-Gaussianities 

The non-commutativity also affects the non-Gaussian 
statistics of primordial perturbations. Here we discuss 

non-Gaussianities in effective single field models, de- 
scribed by (1), where the field (f) affects very little the 
classical dynamics during inflation, V /SH'^M'^ 0, 
but becomes dynamically relevant at some later stage. 
Such models can generate observable non-Gaussianities 
in the usual commutative case (Enqvist & Vaihkonen 
2004; Lyth & Rodriguez 2005), see e.g. works concern- 
ing the curvaton scenario (Bartolo et al. 2004; Enqvist 
& Nurmi 2005) or modulated reheating (Enqvist et al. 
2005; Ichikawa et al. 2008), and our aim is to analyze 
how the non-commutativity alters the predictions. 

Assuming standard slow roll dynamics with canonical 
kinetic terms during inflation, the Fourier transform of 
(1) can expressed in the form 

Ce{t,k) = Ce{t,k) + f{k) j -^le{t,q)C9{t,k-q) + . . . , 

(59) 

where we have defined 



N" 



1 + '^/nl.o In ■ 



N' f ^ V"'\ 



(60) 
(61) 
(62) 



Here we have neglected slow roll corrections to constant 
terms and derived the fc-dependent terms to leading or- 
der precision in slow roll. The same precision is used 
in the results derived below. N = N{t,ti) and all other 
quantities in (60), (61) and (62) without explicit time in- 
dices are evaluated at the time ti appearing in (1). The 
slow roll parameter rj is defined hy r] = MpV" /3H^ and 
V denotes the potential of (p. 
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In the commutative case, the field is Gaussian to 
leading order in slow roll. The magnitude of primordial 
non-Gaussianities is then controlled by the function f{k) 
which is related to the non-linearity parameter /nl for 
equilateral configurations. The scale-dependence of /nl 
is measured by n/^^- 

In non-commutative theories, (0 becomes non- 
Gaussian due to the inherent non-Gaussianities of the 
fluctuations 5(f)e. This affects both the magnitude of 
/nl and its scale-dependence. Quite generally, the non- 
Gaussianities also deviate from the (quasi-)local form 
since the non-Gaussianity in (60) is not of the simple 
(Gaussian)^ type. Below we analyze non-Gaussianities 
arising in non-commutative theories discussed in section 
3. (The non-Gaussianity in the approach of section 4 
can be considered starting from (38), but for clarity we 
do not consider that separately here). 

Using the relations between n-point functions of S(f)o 
and d(po given in section 3, we can express the three- 
point function of (^0 in the form 

{CeitMKeitMCeitM) = i'^^fSiJ^h) 



where Pq = N''^H^/2k^ denotes the power spectrum 
in the commuting case. In this section we denote the 

wavevectors by an ovcrbar, fc = k. 

To obtain real-valued results in coordinate space, we 
apply the prescription introduced in (23) to (63) and 
identify the observable three-point correlator with 

(CeCeCe) = a(CeCeCe)M + «(1 - Q!)(CeCeCe)A 

= {2TTf6{Y,~ki)Be{ki,h,h), (64) 

where the self-adjoint part is 

(Ce(fci)C«(fc2)Ce(fc3))M= (65) 

^((Ce(fcl)Ce(fc2)C^(fc3)) + (Ce(-fcl)Ce(-fc2)C^(-fc3))*), 

and the antiself-adjoint part is 

{Ce{h)Ce{k2)Ce{h))A= (66) 

^((a(fci)Ce(fc2)Ce(fc3))-(Ce(-fci)Ce(-fc2)C«(-fc3))*), 

and all perturbations are evaluated at the same time t. 
It becomes then straightforward to derive the result 

Be{h,k2,h)= (67) 
2e-ifciAfc2 cosh(2if^ • h)Po{ki)Po{k2)f3{h) + 2p.) 

+2ieflf' ( sinh(2//e"*- k3)Po{ki)Poik2)f3{h) + 2p.) . 

We used a shorthand notation +2 p. to denote the per- 
mutations of the three indices. Using the above result 
and the spectrum of the two-point function computed in 
section 3, we find the non-linearity parameter /nl.s given 
by the expression 



/i 



NL,( 



5 Bg(fci,fc2,fc3) 

6Pe{ki)Pe{k2) + 2p. 



= -exp^(-ifci A ^2) 

Po(fci)Po(fc2)/(fc3)cosh(2gg ■ fca) + 2p. 
^ Po{ki)Po{k2)e^iH9 ■ ki)ea{H0 ■ h) + 2p. 
5 

-|-i-expi_^(ifci Afe) 

Po(fci)Po(fc2)/(fc3)sinh(2gg"'. h) + 2p. 
^ Po(A;i)Po(fc2)e„(F^-fci)e„(M-fc2) + 2p. ' 

We have separated the real and imaginary terms in (67) 
and (68). The imaginary contributions violate parity and 
they vanish if all the components of are set to zero. 
Otherwise they arc present for arbitrary a: in particu- 
lar we observe that restriction to the self-adjoint piece of 
the correlation (a = 1) does not eliminate the odd corre- 
lations. One also expects that higher order correlations 
would exhibit parity violations even in the case of purely 
spatial noncommutativity. The spatial components of 
the non-commutativity matrix 9ij enter the through re- 
sults the phase exp(zA:i A fc2) = exp(A:Jfc2 0^ ). They do 
not appear in the results for the spectrum and there- 
fore affect only the non-Gaussian statistics of primordial 
perturbations. 

For simplicity, we analyze in the following only the 
modifications due to Oij setting all components of equal 
to zero. This gives 



/nl,6i — ^e, 



,^^^ Po(fcl)Po(fc2)/(fe3) + 2p. 

Po(fci)Po(fc2) + 2p. 



(69) 



where the only contribution from the non-commutativity 
is the prefactor involving the wedge product. This af- 
fects the scale dependence of /nl,9 and can hence be 
constrained observationally. For example, computing the 
scale-dependence for shape preserying variations of the 
momentum space triangle, fcj Xki, defined as 



51n I /nl,0 ( Afci , Xk2 , Afcs ) I 



ainA 

we find, in the two specific cases, 



(70) 



2k-j^k'20^j cot ( /l -[^ /l2 ) 
2k-^k^Q'ij t an(/L 2 /l2 J ) 



^- "-/nl.q 



if a 
if a 



(68) 



0, 
1. 

(71) 

where w/nl.o given by (62) is the result in the commut- 
ing case. The part dependent on Oij arises purely from 
non-commutative features. The observational prospects 
of scale dependent /nl were considered in Ref. (Sefusatti 
et al. 2009), which suggests that Planck data could be 
sensitive to a scale dependence of the order of slow roll 
parameters. The scale dependence therefore could place 
interesting bounds on 9ij. Moreover, it is worth noting 
that the result (71) depends on the wavevectors ki and 
k2 and hence on the shape of the momentum space tri- 
angle. This is in contrast with the commutative case, 
where the shape dependence is given by the same result 
'^/nlo fo'' all shape preserving variations, ki Xki, re- 
gardless of triangle shape. This allows, in principle, to 
distinguish between the contributions arising from the 
non-commutative properties of the theory and from the 
standard classical inflationary physics. 
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6. CONCLUSIONS 

We considered the effects of noncommutative geometry 
to the statistics of the CMB anisotropy field. The results 

are encoded in two formulas: 

• The statistically anisotropic modulation of the two- 
point function: Eq.(25). 

• The function fj\[L characterizing the non-Gaussian 
property of the three-point function: Eq.(68). 

Both of these describe effects that in general violate 
parity. The presence of spacetime noncommutativity 

was found to induce the leading contributions to the 
anisotropic couplings, which occur in principle between 
all pairs of multipolcs. The non-Gaussianity is scale- 
dependent in a way which depends upon the shape of 
the momentum triangles considered. These features can 
provide stringent bounds on the scale of noncommuta- 
tivity. 

The first tests of these predictions are underway 
(Groeneboom et al. 2010). A quite promising result is 
that already the leading contribution, a dipole modu- 
lation, is found to have an anomalous signature which 



exhibits a hemispherical asymmetry and is modestly pre- 
ferred by the data. We hope to make progress also on 
the theoretical problem of the physical part of the corre- 
lations in a future publication. In particular, the correct 
value of a, introduced in (23), should be deduced from 
first principles, whereas it was left here as an additional 
parameter to be determined empirically. 

By looking closely at the odd features in the sky, one 
may see evidence that an accurate description of the uni- 
verse must be deformed and twisted, since it is fundamen- 
tally pointless. 
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